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Abstract 

The modified model emission function S(x,p), with a uniform radial expansion mode and mass 
reduction, is used to study the single particle spectra and two particle correlation of pions. The 
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calculated HBT radii and one particle spectra are compared with the data from ^s NN = 130 
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GeV Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). We see the main features 
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of the data are reproduced, R /R s decreases with increasing Kt and we have a kinetic freeze out 
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temperature T = 150MeV . 
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I. INTRODUCTION 

We believe quarks are constituents of hadrons, but at present we have not found "free 
quarks" yet. At low densities quarks are confined in individual hadrons, we expected that 
at high energy density nuclear matter undergoes a transition to deconfmed " quark matter" 
or " quark gluon plasma" , where quarks could move " freely" in a volume much larger than 
that of a single hadron. The primary goal of the ultrarelativistic heavy ion collision is to 
create the new form of matter, quark gluon plasma. Bose-Einstein correlation in multipar- 
ticle production processes provide valuable information on the space-time dynamics of the 
fundamental interactions. It is expected that correlations of identical pions produced in the 
ultrarelativistic heavy ion collisions could lead to a better understanding of the properties 
of the new matter form. The main aim of particle interferometric methods is to extract 
as much information as possible about the emission function S(x,K), which characterizes 
the particle emitting source created in the heavy ion collision. A completely model inde- 
pendent reconstruction of emission function from measured correlation data is not possible, 
only certain combinations of spatial and temporal source characteristics are measurable [l| , 
and some useful pieces of information about the phase-space distribution S(x, K) could be 



obtained . One particle spectra anc 
functions are related to it in J2, ua, u \ 

AN 



two particle HBT (Hanbury-Brown-Twiss) correlation 



E lfp = J d ' xS ^P^ W 
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where K — (pi + £>2)/2 ,q = (pi — p 2 ) and p%, p 2 are momentums of detected particles. 
Experimental measurements of C are usually parameterized in terms of the intercept X(K) 
and the HBT radii R%{K) by 

C(K, q) = 1 + A exp [-Y,RU K )ml (3) 

The relative momentum is decomposed into components parallel to the beam 
(l=longitudinal), parallel to the transverse component of K (o=out), and in the remain- 
ing third direction (s=side). A basic model for S(x,K) is proposed in ref. |l|, |a|. We expect 
that its main characteristics can be quantified by its widths in the spatial and temporal 
directions, and a collective dynamical component. 
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Here r = \Jt 2 — z 2 denotes the longitudinal proper time and 77 = |ln[(£ + z)/(t — z)\ the 



space-time rapidity. The transverse mass m± = Jp\ + m 2 . Spin degeneracy (2 J + 1) = 1 
, and chemical potential fi = for pion. The particle four-momentum can be expressed by 
the momentum rapidity and the transverse mass, 

Pfi = {m ± coshy,p ± ,0,m ± smhy), (5) 

p-u(x) m± p±x . 

T =—cosh(y-r])coshr] t -—-smhr] t . (6) 

The source has a finite geometrical size in the spatial and temporal directions encoded in 
transverse and longitudinal Gaussian widths R and A?? as well as in a finite particle emission 
duration At. The freeze out temperature are mainly determined by the the transverse single 
pion spectra p and resonance decay contributions are often taken into account in fitting 
processes 0] . The influence of collective velocity is considered also with a linear transverse 



flow rapidity profile 



Vt(r) = rif— (7) 



in the model [8J]. With these assumptions the HBT radii can be calculated by 

R 2 (K) =< f > (K), (8) 

R 2 (K)=<(i-(3j) 2 >(K), (9) 

R 2 (K)=<(~z-p l i) 2 >(K) (10) 

where < ... > denotes an average with the emission function S(x,K): 

<^>w= j ffi^r ) - (n) 

The space-time coordinates x M are defined relative to the "effective source center" x^ by 

x^{K) = x^-x^{K), 

x^(K) =< x^ > (K). (12) 



But at present it is difficult to reproduce one particle spectra and the HBT radii extracted 
from two-pion correlations simultaneously. R /R s from hydrodynamics is greater than one, 
and increases as K? does. In stead, experiment finds that R / R s decreases with increasing 
Kt [9] . Different approaches are proposed to deal with this problem recently |l£], |ll|, |l2| . 
Here we would like to modify this basic model and try to find a set of parameters, by which 
the main features of the one particle spectrum and HBT radii could be reproduced. 

II. THE MODIFIED MODEL 

(l)We use Bose-Einstein distribution function instead of Boltzmann function. The differ- 
ence is not small for pion at low transverse momentum. The resonance decay contributions 
are not included. Only the direct pions are considered, 
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(2)The radial expansion mode is assumed to be uniform, independent of radial position, 
i] t = C. We are motivated to do this by the fact that the experiment shows a type of 
"explosive" behavior, which implies that a strong pressure is build-up at the origin. We first 
generalize the linear transverse flow rapidity profile to 

Vt(r) = Vf(^) a - (14) 

It is easy to see, a smaller a implies a stronger radial flow. Then we take a = to express 
the strongest radial flow, the strongest pressure build-up at the origin. Maybe a 5 function 
could be used to describe the pressure at the origin. We neglect the detailed structure in 
this small region near the center and express it in a simplified way. 

(3)In the procedures of fitting, the temperature influence on the pion mass is considered. 
Many different models are used to approach this problem, their results are quite different. 
We only use the effective pion mass as an adjustable parameter here, its value is fitted by 
the data. We will compare it with some model estimated values in section III. 

Under these modifications, a set of parameters are found, m^ = 40MeV, T = 150MeV, 
r\ t = 0.55, R = 5.2/m, Ai] = 5.6, r = 6/m/c, and At = lfm/c. 



The Pt dependence of the yield of 7r is fitted by 
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Afd 4 xS{x,P) (15) 



P T dP T dy 

It is assumed that all particles decouple kinematically at the same freeze out temperature. 
The experimental data are taken from ref. [l3|. Transverse momentum spectra of n~ are 
measured at midrapidity for the most central events. The results from the emission function 
with mass reduction, uniform radial expansion and Bose-Einstein distribution are shown in 
fig. 1 by the solid line. It is fitted in a broad range Pt < 2.2GeV. The dashed line shows 
the results from the emission function with the uniform radial expansion and Bose-Einstein 
distribution, but without mass reduction. We see that the influence of mass reduction exists 
in low transverse momentum region. The dot dashed line shows the results from the emission 
function with linear radial expansion and without mass reduction (other parameters are the 
same as given above) . The curve is changed greatly. This implies that the extracted 
temperature is connected with the expansion mode tightly. 

In fig. 2 solid line shows the calculated results of R s by the emission function with mass 
reduction and uniform radial expansion. Data of R s for pions as a function of K? are 
measured at mid-rapidity from y/s NN =130 GeV Au+Au collisions at RHIC. Solid circle 
(empty circle) for tc + (n~) are given by ref. [14J. Triangle-down (triangle-up) for n + (n~) 
are given by ref. 15j. We see the main features of the data are reproduced by the modified 
emission function. The dot dashed line shows the results from the emission function without 
modifications. 

The experimental data for R and Ri are also compared with the calculated results in 
fig. 3 and fig. 4 respectively. The results from the emission function with mass reduction and 
uniform radial expansion are shown by the solid lines and the results from the emission 
function without modifications are shown by the dot dashed lines. We see the slope of R 
is larger in the modified model. The ratio R /R s obtained from the emission functions 
is shown in fig. 5 by the solid line (with modifications) and the dot dashed line (without 
modifications). Data are from PHENIX [14j and STAR [l5j. The pion source parameters 
given here do not reproduce the data exactly yet, but it is consistent with the main features 
of the single particle spectra and HBT radii, R /R s decreases with increasing K?- If this 
picture is acceptable, we have a kinetic freeze out temperature T = 150MeV. It is just near 
the chemical freeze out temperature. 



III. DISCUSSIONS 

(l)We take the uniform radial expansion mode here. It does not mean that there is only 
pressure (and pressure gradients) at the origin of the system. It only means that pressure 
gradients at the small region near the origin of the system is much stronger than that at 
the other parts. It is only a simplified description. We take this choice as it can produce a 
reasonable result and it is consistent with the picture of explosion. 

(2)The kinetic freeze out temperature is mainly determined by the one particle spectra. 
With different expansion mode the extracted temperature is different. In fig.l the dashed line 
shows the result from the emission function with linear expansion profile and without mass 
reduction. It is only used to show the influence of the modification. Our result T = 150MeV 
is just near the chemical freeze out temperature. It implies that the chemical freeze out and 
kinetic freeze out maybe occur simultaneously. 

(3)A basic idea for deconfinement is that the vacuum can have two possible states. The 
first is the normal vacuum outside the hadrons. Quarks and their gluon fields modify the 
vacuum in their neighborhood, transforming it into the second state, the " perturbative" 
vacuum, or "bag" state, the vacuum inside the hadron. The inside state of the vacuum 
has a higher energy per unit volume than the outside vacuum state, the energy difference 
, B, is often called as bag constant in the bag model. The current quark mass is about 10 
MeV for u quark and d quark. The mass of hadron with light quarks, i.e. the energy of a 
static spherical bag with radius R, is E = ^ + BR where V is the volume of the bag and 
C is a positive constant, which is determined by the number of quarks in the bag and the 
eigenvalues of energy and zero-point energy [la,ll2|. Determined by dE/dR = 0, the energy 
of the bag has a minimum value E = 4/3(4-7nB] C 1//4 . R.D.Pisarski gives the temperature 
dependent bag constant _B(T)oc(T 4 — T 4 ) in [18[. The field theory at finite temperature also 



shows the bag constant is temperature dependent [19|, |20|, B(T) = B (l — |^) 2 where T c 
is the critical temperature. At the critical temperature the difference between the normal 
vacuum and the perturbative vacuum disappears, and the quarks are deconfinement. Near 
the critical temperature the difference between the two vacuums is small and thus the mass 
of hadron is reduced. At the freeze out condition, the hadrons decouple from each other, 
the mass of a hadron really has its vacuum value, but as the temperature around it is still 
high, it has a "hot vacuum" value. A "free hadron" in the vacuum is not "alone", the 



vacuum around it is full of virtual particles. For hadrons with heavy quarks, H. Satz use 



potential at finite temperature to study their properties 



21| . The detailed expressions for 



mass shift are different in these approaches, but they have common points. The hadron 



mass decreases with the increasing temperature 



The mass shift is also studied in other 



ways. R. D. Pisarski and F. Wilczek |23( study the restoration of the spontaneously broken 
chiral symmetry in a linear a model and show that pion mass is stable but under some 
conditions the mass of 7i"o could drop to 77 MeV. Pion is taken as Goldstone bosen its mass 
would be zero, but in the realistic case, pion mass has a small value. G. E. Brown and M. 
Rho studied the chiral restoration in hot and/or dense matter, and showed that the explicit 
chiral symmetry breaking-that gives masses to the Goldstone bosons could be "rotated 
away" in J24I] . We also use linear a model to study the properties of 7r° J25|, we find that 
if the explicit breaking term is temperature independent pion mass could even increase as 
temperature does, but the temperature dependent coupling will change the situation. From 
the above discussions we see the predictions from different models are quite different. O nly 
the experiments could give judges. Unfortunately we have not found enough evidence [2f| 
to show that the hadron mass (and pion mass) decreases with increasing temperature at 
present. It is difficult to measure the hot hadron mass directly as the observer is in low 
temperature environment, the pion mass returns to the normal value there. From the view 
point of bag model, the volume of pion shrinks and obtains the energy from the vacuum 
(bag constant), thus the pion mass returns to the normal value. In this process, the inter 
structure of pion is changed. Such a change will give no influence on the Pt distribution of 
pion. We could measure the influence of mass shift , but in some cases the influence could 
be even washed away by other facts. In our fitting procedures, we only take pion mass as an 
adjustable parameter to fit the data. Our result seems consistent with the mass reduction 



picture and previous studies |27, 



28, 



29 



3l|. 



(4) The correlation function is determined by the momentum distributin at the freeze 
out temperature in both the original model and the modified model. The difference of the 
momentum distributions at the freeze out temperature between these two models are caued 
different thermal distributions before the freeze out. 

(5)A completely model-independent reconstruction of the emission function S(x,p) is not 
possible. We follow the approach given in refs. [l|, |5| by using simple parametrizations of the 
emission function. We adjust the parameters by a comparison to data. The modified model 



can reproduce the main features of the one particle spectra and HBT radii simultaneously 
now. This model should be further improved, if it is used to describe the azimuthal mo- 
mentum space anisotropy of particle emission [32[. The pion production from resonance is 
not considered here. In ref.(l|, resonance decay contributions are shown. Thus a detailed 
calculation including these effects is necessary for us to make the model better. We will 
discuss it elsewhere. 
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Captions: 

Fig.l Transverse momentum spectra measured at mid-rapidity for 7r~ at most centrality 
given in ref.yjj] are compared with results obtained from the emission function with mass 
reduction, uniform radial expansion and Bose-Einstein distribution(solid line), with uniform 
radial expansion and Bose-Einstein distribution, but without mass reduction (dashed line), 
and with linear radial expansion and Bose-Einstein distribution but without mass reduction 
(dot dashed line). 

Fig. 2 R s for pions as a function of K T are measured at mid-rapidity. Solid circle (empty 
circle) for 7r + (7r~) are given by ref. [14|. Triangle-down (triangle-up) for n + (n~) are given 
by ref. 15(. Solid line shows the calculated results from the emission function with mass 
reduction and uniform radial expansion. Dot dashed line shows the results from the emission 
function without modifications. 

Fig. 3 R for pions as a function of K T are measured at mid-rapidity Solid circle (empty 
circle) for n + (n~) are given by ref. [14]. Triangle-down (triangle-up) for n + (n~) are given by 
ref.|l^|. Solid line shows the calculated result from the emission function with mass reduction 
and uniform radial expansion. Dot dashed line shows the results from the emission function 
without modifications. 

Fig. 4 Ri for pions as a function of Kt are measured at mid-rapidity Solid circle (empty 
circle) for n + (n~) are given by ref. |l4j. Triangle-down (triangle-up) for n + (n~) are given 
by ref. |l£| . Solid line shows the calculated result from emission function with mass reduction 
and uniform radial expansion. Dot dashed line shows the results from the emission function 
without modifications. 

Fig. 5 Rq/R s for pions as a function of K T are measured at mid-rapidity Solid circle 



(empty circle]_for ix + (ix ) are given by ref. 14]. Triangle-down (triangle-up) for tt + (jt ) are 



given by ref.(l5j]. Solid line shows the calculated result from emission function with mass 
reduction and uniform radial expansion. Dot dashed line shows the results from the emission 
function without modifications. 
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